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The response of a quantum system to a static field has 
been a longstanding problem since the early days of quan- 
tum mechanics. A topic of particular interest in this wide 
field is the dynamics of a quantum particle in a periodic 
potential induced by a static force (modelling a crystal 
electron in an electric field). In this system, the effect of 
the field manifests in a very unintuitive way Indeed, as 
already emphasised by Bloch Tj and Zener |2| , according 
to the predictions of wave mechanics, the motion of elec- 
trons in a perfect crystal should be oscillatory rather than 
uniform. This phenomenon, nowadays known as Bloch 
oscillations (BO), has recently received renewed interest 
which was stimulated by experiments on cold atoms in 
optical lattices 0, S IE 13 • This system (which mimics 
a solid state system - with the electrons and the crystal 
lattice substituted by the neutral atoms and the opti- 
cal potential, respectively) offers unique possibilities for 
the experimental study of BO and of related phenom- 
ena. In turn, these fundamentally new experiments have 
stimulated considerable progress in theory (see review 
0, and references therein), and it can be safely stated 
that BO in diluted quasi one-dimensional gases is well 
understood today. Other directions of research focus on 
BO in the presence of relaxation processes (spontaneous 
emission) |g, BO in 2D optical lattices 0, and BO in 
the presence of atom-atom interactions ('BEC-regime') 
|i0,IL jJ,lll- The present Letter deals with the third 
problem, which is approached here by an 'ab initio' anal- 
ysis of the dynamics of a system of many atoms. This 
distinguishes this work from previous studies of BO in 
the BEG regime 0, 0, , which were based on the a 
mean field approach using a nonlinear Schrodinger equa- 
tion. A new effect, so far unaddressed by these earlier 
studies, is predicted: besides the usual Bloch dynamics, 
the atomic oscillations may exhibit another fundamental 
period, entirely defined by the strength of the atom-atom 
interactions. 

Let us first recall some results on BO in the single- 
particle case. Using the tight-binding approximation (14 | , 
the Hamiltonian of a single atom in an optical lattice has 
the form 
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In Eq. \l) denotes the ^th Wannier state (/);(x) cor- 
responding to the energy level Eq 0, J is the hopping 
matrix elements between neighbouring Wannier states, 
d is the lattice period, and F is the magnitude of the 
static force. The Hamiltonian can be easily diag- 
onalised, which yields the spectrum Ei — Eq -\- dFl 
(the so-called Wannier-Stark ladder) and the eigenstates 
(Wannier-Stark states) 

\A) = ^ J™-/(J/dF)|m) , {x\m) = 0„(x) , (2) 



(here Jm{z) are the Bessel functions). As a direct con- 
sequence of the equidistant spectrum, the evolution of 
an arbitrary initial wave function is periodic in time, 
with the Bloch period Tb = 2Trh/dF. In particular, 
we shall be interested in the time evolution of the Bloch 
states I ■0k) = X]/ exp(zdKZ)|Z). Using the exphcite expres- 
sion for the Wannier-Stark states l(2Jl, it is easy to show 
that Itpnit)) = exp{— i( J/dF) sin(dK(t))}|?/'K(t)), where 
K{t) = K -\- Ft/h (from now on iSg = for simplicity). 
Note that the exponential pre-factor in the last equation 
contains the same parameter J/dF as the argument of 
the Bessel function in Eq. Depending on the value 
of this parameter, the regimes of weak {dF <C J) or 
strong {dF ^ J) static fields can be distinguished. In 
this Letter, we shall restrict ourselves to the strong field 
case, which, in some sense, is easier to treat than the 
weak field regime. Indeed, for J/dF <C 1, the Wannier- 
Stark states practically coincide with Wannier states, and 

\Mt)) « IV-^w). 

A remark concerning the characteristic values of the 
parameters is at place here: In the numerical simulations 
below, we use scaled variables, where ?i = 1, d = 27r, 
and the energy is measured in units of the photon re- 
coil energy. In typical experiments with cold atoms in 
an optical lattice, the amplitude v of the optical poten- 
tial equals few recoil energies. Then, for example, for 
V equal to 10 recoil energies, the value of the dimen- 
sionless hopping matrix element is J = 0.0384. The 
strength of the static field is restricted from below by 
the condition dF > J, and from above by the condition 
that Landau-Zener tunnelling events can be neglected. 
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FIG. 1: Momentum distribution of the atoms in the opti- 
cal lattice, for different amplitudes v of the optical potential. 
(The amplitude v is measured in units of the recoil energy, 
the momentum k in units of 2-Kfi/d.) The figure illustrates 
the transition from the SF-phase to the Ml-phase as v is var- 
ied (F = 0, L = = 7). 



Since the probability of Landau-Zener tunnelling is pro- 
portional to exp(— 7r(5^/8(iFJ) [5 is the energy gap sepa- 
rating the lowest Bloch band from the remaining part of 
the spectrum) 0,13, we have F < 30 for w = 10. 

We proceed with the multi-particle case. A natural 
extension of the tight-binding model , which accounts 
for the repulsive interaction of the atoms, is given by the 
Bose-Hubbard model jl^ , 

\/=i / 1=1 
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In Eq. |(2Jl, and di are the bosonic creation and anni- 
hilation operators, fii = ajai is the occupation number 
operator of the Ith lattice site, and the parameter W 
is proportional to the integral over the Wannier func- 
tion raised to the fourth power. Since the Bose-Hubbard 
Hamiltonian conserves the total number of atoms N, the 
wave function of the system can be represented in the 
form 1 4') ~ J2ii'^n\n), where the vector n, consisting of 
L integer numbers ni n; = N), labels the iV-particle 
bosonic wave function constructed from N Wannier func- 
tions. (In what follows, if not stated otherwise, \'^) refers 
to the ground state of the system.) As known, in the ther- 
modynamic limit, and for F = 0, the system ^ shows 
a quantum phase transition from a superfluid (SF) to a 
Mott insulator (MI) phase as the ratio J/W is varied (see 
(l?! and references therein). It is interesting to note that 
an indication of this transition can already be observed 



FIG. 2: Bloch oscillations of the atoms, induced by the static 
force F = l/27r (v = 10). One Bloch period is shown. 

in a system of few atoms |0. As an example. Fig. 
shows the diagonal elements of the one-particle density 

matrix, 

L 

p{k,k') = (*|$t(fc)$(fc')|^) , ^k) ^Y^^aMk) , 

1=1 
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for TV = L = 7, 5 < i; < 35, and = 0.1 / dk(t>j{k) 
(here, 4>i{k) are the Wannier states in the momentum 
representation, and k — p/{2T:h/d) is the dimcnsionless 
momentum). Physically, this quantity corresponds to the 
momentum distribution P{k) — p{k, k) of the atoms, di- 
rectly measured in the experiment. It is seen in Fig. ^ 
that, around v = 15, there is a qualitative change in 
the momentum distributiori, in close analogy with that 
observed in the experiment [13 • It should be noted, how- 
ever, that this qualitative change of the momentum dis- 
tribution alone does not yet prove the occurrence of a 
phase transition. A more reliable indication of a SF-MI 
phase transition are the fluctuations of the number of 
atoms in a single well, which drops from (An^) w 0.72 at 
w = 5 to (An2) « at u = 35 

Let us now discuss the effect of the static force. Figure 
121 shows the dynamics of the momentum distribution of 
the atoms (which were initially in SF-phase) in presence 
of a force F — 1/2tt This numerical simulation illus- 
trates atomic BO as observed in laboratory experiments 
0, 01 . It is seen that after one Bloch period the initial 
momentum distribution practically coincides with the fi- 
nal distribution. A small difference, which can be noticed 
by closer inspection of the figure, is obviously due to the 
atom-atom interaction |2^ . This difference becomes evi- 
dent once the system evolved over several Bloch periods. 
In Fig. 121 the momentum distribution P{k) at integer 
multiples of the Bloch period is shown. A periodic change 
of the distribution from SF to Ml-like and back is clearly 
seen. (The use of the term 'Ml-like' stresses the fact that 
the variance (An^) does not change as time evolves.) In 
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FIG. 3: Dephasing of Bloch oscillations due to the atom- 
atom interaction. The period Tw = 2nF/W is clearly seen. 
(F = l/2n, v = 10,W = 0.1 / dx(j)f{x) = 0.0324.) 
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FIG. 4: Normalised mean momentum (p/NJ — » p) as a func- 
tion of time, for two different values of the occupation number 
n = N/L: N = L = 7 (upper panel), and iV = 4, L = 14 
(lower panel) . The dashed lines show the analytical result for 
the envelope function in the thermodynamic limit. 



addition to Fig. |21 and Fig. |21 Fig. ^ depicts the mean 
momentum p{t) of the atoms for two different values of 
the occupation number n = N/L (number of atoms per 
lattice cite) - n ~ 1 (upper panel) and n — 2/7 (lower 
panel). As to be expected, the dynamics of the system 
depends on the value of n, and for a larger occupation 
number the deviations of many-particle BO from the non- 
interacting result p{t) — NJ sm{2TTFt) becomes larger. 

Our explanation for the numerical results is the follow- 
ing. It is convenient to treat the atom-atom interaction 
as a perturbation. Let us denote by Uo{t) the evolution 
operator of the system for VF = 0, by U{t) the evolution 
operator for W ^ 0, and by Uwit) the evolution operator 
defined by the decomposition U{t) = Uw{t)Uo{t). Since 



Uq{Tb) is the identity matrix, one has to find UwiTB) to 
reproduce the result of Fig. |31 In the interaction repre- 
sentation, the formal solution for UwiTB) has the form 



Uw{Tb) = exp 
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where the hat over the exponential denotes time order- 
ing. Now we make use of the above strong-field condition 
dF > J. Under this premise the Wannier states are the 
eigenstates of the atom in the optical lattice [see Eq. 0]. 
In the multi-particle case this means that the Fock states 
|n) are the eigenstates of the system (jSJ and, thus, that 
the operator Uo{t) is diagonal in the Fock state basis. 
Then the integral in Eq. jSJ can be calculated explicitely, 
which yields 
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{n\Uw{TB)\n) = exp }_^{n\hi{ni - l)|n) j . 

(6) 

Finally, by noting that the quantity {n\ni{ni — l)|n) is 
always an even integer, one comes to the conclusion that, 
besides the Bloch period, there is additional period. 



Tw = 2ttF/W 
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which characterises the dynamics of the system. 

Further analytical results can be obtained if we approx- 
imate the ground state of the system for = by the 
product of N Bloch waves with quasimomentum k = 0, 
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Indeed, let us consider, for example, the dynamics of the 
mean momentum. Using the interaction representation 
(now with respect to the Stark energy term) the mean 
momentum is given by 

Pit) = Jim (^Wlit)\Y^^a\^,di\Uwim)e-'^^^'^ , 

where Uwit) is the continuous- time version (Tg — )■ t) 
of the diagonal unitary matrix © . Substituting Eq. ^ 
and Eq. 10 into Eq. 0, we obtain the following exact 
expression. 
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where Vin,n') is the joint probability to find n and 
n' atoms in two neighbouring wells. In the thermo- 
dynamic limit N,L — )■ oo, N/L = n, the function 
Vin^n') factorises into a product of the Poisson distri- 
butions Vin) = n" exp(— n)/n!, and the double sum 
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in Eq. (|10|l converges to the positive periodic function, 
f{t) = exp(-2n[l - cos{Wt)]), indicated in Fig. H by 
the dashed hne. Good agreement between the envelope 
of p{t) and the dashed hne proves that in the numerical 
simulation presented above the convergence was indeed 
achieved. 

In summary, Bloch oscillations of interacting cold 
atoms have been studied, both numerically and analyt- 
ically. We have shown that in the strong field regime 
atom-atom interactions cause the reversible dephasing 
of Bloch oscillations. As a result, the momentum dis- 
tribution of the atoms changes periodically from SF to 
Ml-like distributions, with a period given by Eq. 10 . Us- 
ing the original (unsealed) parameters, this period reads 



Tw = {dF/W)TB = 21:%/ W, where W is the strength 
of the atom-atom interactions. Since the momentum 
distribution can be measured easily in the laboratory 
experiment, this effect suggests an alternative method 
for studying atom-atom interactions by applying a static 
force to the system. It is worth to stress one more time 
that the reported result is valid only in the strong field 
limit, dF > J. If this condition is violated, the evolu- 
tion operator © is no more a diagonal matrix, and the 
system dynamics get significantly more complicated. An 
analysis of this latter case will be presented elsewhere. 

Discussions with A. Buchleitner are gratefully acknowl- 
edged. 
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